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Abstract We report on a lattice investigation of improved quenched Wilson fermions 
above and below the confinement-deconfinement phase transition. Results on meson 
screening masses as well as spatial wave functions are presented. Moreover, the meson 
dispersion relation is studied. Below the critical temperature we do not observe any 
significant temperature effect while above T c the data are consistent with a leading 
free quark behavior. 



1 Introduction 



An important goal of analytical as well as lattice investigations of QCD at non- 
vanishing temperature has been to gain more insight into the temperature depen- 
dence of hadron properties below and into the nature of hadronic excitations above 
the transition temperature from the hadronic to the plasma phase of QCD. 

When the temperature is raised towards the transition point, approaching (ap- 
proximate) chiral symmetry restoration and deconfinement is expected to change the 
properties of hadrons. In particular, the lightest vector mesons and the temperature 
dependence of their masses and decay widths have received quite some attention be- 
cause of their possible relevance to the observed enhancement of low mass dileptons in 
heavy ion collisions 0. The theoretical predictions of these properties are, however, 
model dependent, see e.g. for a recent review. In the plasma phase, the effective, 
temperature dependent coupling constant becomes small at large temperatures. One 
is thus lead to expect that the plasma consists of a gas of only weakly interacting 
quarks and gluons. In this case, correlation functions of operators with hadron quan- 
tum numbers should be described by the exchange of two or three almost free quarks. 
On the other hand, there are arguments that even at high temperature the hadronic 
excitation spectrum might be more complicated because of non-perturbative effects 
in particular in the chromo-magnetic sector of QCD ||. Thus, in both temperature 
regimes ab-initio QCD computations are highly desirable. 

Hadronic correlation functions at non- vanishing temperature T have been the sub- 
ject of lattice investigations for quite a while. Most of these studies are based on the 
staggered quark formulation @-14] while only few so far utilized the Wilson discretiza- 



tion [p~5|— pr7[1 and, just recently, domain wall fermions |18|]. Moreover, because of the 
limited extent of the lattice in the Euclidean time direction, < t < 1/T, most of 
these studies investigated spatial correlation functions and screening masses. 

Below the transition temperature, detailed lattice investigations and comparisons 
of hadron masses at T = and at T < T c have been carried out in the quenched 
approximation so far. Analyses in the staggered fermion discretization [12[| have 
covered temperatures between 0.8T C and 0.95T C and did not detect any significant 
difference between zero-temperature masses and finite-temperature screening masses 
in the investigated hadron channels. This is also supported by a recent study of spatial 
and temporal correlation functions of Wilson fermions on anisotropic lattices [|17]] at 
T ~ 0.93T C . 

At temperatures above T c , the available lattice results on hadronic correlation 



functions reflect chiral symmetry restoration because the masses [|17]] and screening 
masses IfSHSj |10|. |13|, |l4f in the vector (p) and axial vector (ai) channel become degen- 
erate. In the pseudoscalar and scalar channels it is observed that it and fo/a become 



degenerate at T c |TJj while the mass difference between n and a /5 seems to van- 
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ish only at higher temperatures, see in particular [|T^] for a chiral extrapolation at 
temperatures closely above T c . These findings are in accord with the expected restora- 
tion of SUr(N f ) x SUl(N f ) and indicate that the anomalous Ua(1) symmetry is not 
effectively restored at the critical temperature. 

As far as the values are concerned, in the staggered discretization vector and axial 
vector screening masses are compatible with the prediction of lowest order perturbation 
theory i.e. the propagation of (almost) free quarks. Scalar and pseudoscalar channels, 
however, show substantial deviations from this expectation, at least in the temper- 
ature interval between T c and 2T C . One might argue that this observation indicates 
the existence of (pseudo) scalar bound states. However, the lowest mass hadron is 
unlikely to change from a mesonic state to a quark-like quasi-particle at a non-critical 
temperature. In addition, studies which apply different boundary conditions [|TTJ sug- 



gest that the spectral function is dominated by a two quark cut. For Wilson quarks 



it is observed |17j that already at a temperature of about 1.5T C the pseudoscalar and 
vector screening masses are very close to each other. There it also is found that a near 
degeneracy holds true for the masses. Depending on the source operator utilized, the 
pseudoscalar is sometimes even heavier than the vector meson. While this behavior 



can at least qualitatively be explained by the propagation of (almost) free quarks, fT7 
provides conflicting evidence since their study of the wave functions, on the other 
hand, indicates the presence of bound states. 

In addition to these observables, also spatial wave functions have been analyzed 
|19| . Here one has found a similar behavior as at zero temperature. The observed 



exponential decay has then been taken as to suggest that the relevant hadronic excita- 
tions are bound states also in the plasma phase, at least at temperatures just above T c . 
According to |21J , this behavior could, however, also be explained by the fact that the 
dimensionally reduced, 3-D effective theory and correspondingly spatial Wilson loops 
in 3+1 dimensions show confinement [plH 23"H- Solving a two-dimensional Schrodinger 
equation with a potential which includes a temperature dependent (spatial) string 
tension leads to exponentially decreasing spatial wavefunctions. The corresponding 
effect on the screening masses [El] would be an 0(10%) correction at the investigated 



temperatures [23[ which so far could not yet be checked quantitatively. 



None of the lattice investigations of hadronic masses at non- vanishing temperatures 
has attempted to carry out the continuum limit. Most of the mentioned analyses are 
based on the staggered discretization. A straightforward computation in the Wilson 
formulation of lattice QCD and a comparison of the results with the staggered ones 
would thus help to gain an idea about the discretization effects. This is the main goal 
of this paper. 

Since we are extracting screening masses and Lorentz invariance is lost at finite T 
due to the heat bath, we have also computed spatial correlation functions projected 
onto some non- vanishing momenta and on the lowest non- vanishing bosonic Matsubara 
frequency. The purpose of this attempt is to test for a sizeable difference between 



3 



spatial and temporal momentum contributions to the vacuum polarization tensor. 

As a by-product of the attempt to construct meson operators with good overlap 
to the groundstate we also were able to obtain information about the Bethe-Salpeter 
amplitudes of the investigated mesons, the pion and the p. 

The paper is organized as follows: in the next section we present some details of the 
simulation. This is followed by the presentation and discussion of the obtained screen- 
ing masses in the pion and the p channel, both below and above the deconfmement 
transition. In section |] we test the dispersion relation at non-vanishing temperature. 
Section ^| contains our estimates of the wave functions. 



2 The simulation 

The results to be presented here are based on gauge field configurations which have 
been generated with the standard Wilson gluon action. We used a pseudo-heatbath 
algorithm [p6j with FHKP updating [27] in the SU(2) subgroups. Each heatbath it- 



eration is supplemented by 4 overrelaxation steps |28|. We have simulated at three 
values of the bare coupling, (3 = 6/g 2 = 6.0, 6.2 and 6.4. At these (3 values the lattice 
spacing has been determined from quite a variety of observables. Depending on the 
quantity, the results spread over a range of about 10 % of the central value, however, 
within the error bars, agreement is observed. In order to obtain the physical temper- 
ature of the lattices in units of the critical one, for definiteness we have consistently 
chosen to set the scale by the string tension f29fl , T/T c = (T/y/a) ■ (y/a/T c ), taking 
T c j yfa from |3(| . The temporal extent of the finite temperature lattices of N T = 8 
then corresponds to physical temperatures of T = 0.93(1)T C , 1.23(1)T C and 1.63(2)T C 
at (3 — 6.0, 6.2 and 6.4 respectively. The spatial volume of the lattices was chosen 
as 24 3 and 32 3 at (3 = 6.0, 24 3 at (3 = 6.2 and 24 2 x 64 at (3 = 6.4. In addition to 
the finite temperature simulations we have carried out runs at zero temperature on 
lattices of size 16 3 x 32 at (3 = 6.0 and 24 3 x 48 at 6.2 mainly in order to supplement 
the available literature data on meson masses by results at quark masses in the range 
between the strange and the charm quark mass. We have generated between 50,000 
and 80,000 gauge field configurations and have analyzed configurations separated by 
500 or 2000 sweeps (see Table [I]). Autocorrelations have been checked to be negligible. 

For the fermion part of the action we used the 0(a) Symanzik-improved Sheikholes- 
lami-Wohlert action E5| with a tree level clover coefficient of c sw = 1.0. The inversion 
of the Dirac matrix was carried out by means of an overrelaxed MR algorithm |3l[ 



at moderate to large quark mass values and by the BiCGStabI for light quarks. 
In both cases an even-odd partitioning was employed. We encountered only very few 
moderately exceptional configurations. Including or omitting those did not modify 
the expectation values. With the same combination of actions as in the present study, 
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the (latest) values of the critical hopping parameter at zero temperature defined by 
the vanishing of the pion mass have been determined as k c = 0.14556(6) |33j and 
k c = 0.14549(2) at (3 = 6.0, k c = 0.14315(2) J3§ and k c = 0.14315(1) 
at (3 = 6.2 as well as k c = 0.14143(3) |34[ at /3 — 6.4. These numbers may be 
used to convert the various values of the hopping parameter into an estimate of the 
corresponding bare quark masses by means of 

This has been used also in Table [l] where our run parameters are summarized. 



Simulation Parameters 




1 d\"\"\r i c^ I 1 I cfcif lcfipc 
J-jdLLlCt; 1 J 1 c bLcxLlbLlCo 


rh fllq [-Lvlt; V J 


6.0 0.105 


16 3 x 32 0(100) 
24 3 x 8 0.93 0(100) 
32 3 x 8 0.93 0(40) 


0.141 200 
0.130 750 
0.128 800 
0.141 200 

U.loU ( OU 

0.128 800 
0.145 30 
0.1445 50 
0.144 80 
0.141 200 


6.2 0.077 


24 3 x 48 0(100) 
24 3 x 8 1.23 0(30) 


0.141 150 
0.130 800 
0.128 950 
0.1428 20 
0.14232 60 
0.14151 100 
0.136 450 
0.128 950 


6.4 0.058 


24 2 x 64 x 8 1.63 0(100) 


0.1409 50 
0.1406 80 
0.1403 110 
0.14 140 
0.13 1100 



Table 1: A summary of our run parameters. The physical value of the quark masses 
has been estimated by using eq. ([[]) together with the physical value of the lattice 
spacing as obtained from string tension determinations. For definiteness, a/cx has been 
taken as 420 MeV. 

From the computed quark propagators we constructed correlation functions of op- 
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erators with the quantum numbers of the pseudoscalar and the vector meson. To 
improve the projection onto the lowest energy state, a gauge invariant extended oper- 
ator was used on the sink site l3f 



M R (x) = ^* i a (x)U i >>(x^x + RZ)r atf ,* j fj (x + RZ) (2) 

±e 

Here, $ and \I> are a quark and an antiquark field separated by a distance R, i,j 
denote color indices and a, (3 are spinor ones. Both indices are to be summed over. 
The quantum numbers are selected by choosing T = 75 for the pseudoscalar and T = 7^ 
for the vector channel respectively. In the vector channel we have averaged over the 
polarizations /i perpendicular to the correlation direction. The explicit sum in eq. @) 
is over all unit vectors e perpendicular to the correlation direction. This extended 
operator is made gauge invariant by introducing the color parallel transporter U from 
x to x + Re. In order to further improve the projection and in an attempt to resemble 
the gluon cloud |36| the parallel transporter is built from smeared link fields. At the 



source site a strictly local operator (R = 0) was put on the lattice. 
The most general correlation function Cn{p,t) is thus obtained as 

c R (p,t) = J2 e ~ ip£ ( p R&t)^&M 

X 

= <0\P R \P(p)><0\P\P(p)>* {e- Ep ^ t + e- Ep ^ (N ^ t) } + --- (3) 

where in this particular example the pseudoscalar correlation in the temporal direction 
was chosen. In eq. (0), the exponential fall-off is given by the energy Ep(p) of the state 
\P{p) > at momentum p. The dots indicate contributions from excited states with the 
right quantum numbers. At non- vanishing temperature, because of the limited extent 
of the lattices in the temporal direction, t < l/T, we computed spatial correlation 
functions in the z direction, Cp(p,z), where p denote the momentum components 
perpendicular to the z direction, p = (p x ,P y ,Pt)- 

The improvement procedure described above leaves quite some freedom in choosing 
optimal parameters. As for the distance R between quark and antiquark, on each 
configuration we have computed the correlation functions for a variety of different R 
values. This allowed us to find the optimal separation for each lattice spacing and 
temperature individually. An example of how the separation R changes the projection 
to the lowest mass state is shown in Figure [|. Here we plot the effective mass, defined 
as M eS (t) = \n{Cji(t)/Cji(t + 1)} for p = 0, of the pseudoscalar as a function of t 
for various R values. The data has been obtained on the zero temperature lattice 
at (3 = 6.2 at k = 0.141. The plot illustrates that the contribution of excited states 
becomes considerably smaller when R is raised from 1 to 5 lattice spacings in this 
example. The data flattens off and reaches a plateau at smaller time separations 
between source and sink. This allows to extract the lowest mass much more reliably. 
At R > 7 the large t limit is approached from below as a single term contributing to 
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eq. fl3|) is not positive-definite. At the rightmost data point in the figure the effective 
mass drops slightly because the periodicity of the lattice is being felt by the correlation 
function. When f3 is varied we observe that the optimal distance in lattice units 
approximately scales with the lattice spacing i.e. stays constant in physical units. 

Contrary to the strong R dependence of the effective masses which appears to be 
physical, details of the gauge field smearing procedure do not seem to matter so much. 
We have adopted the APE prescription |37j with a weight of 2 for the link term and 
1 for the contribution of the staples. As noted also in []36[| , the precise value of the 
ratio is not too important, contrary to the case of Wilson loops. The smeared links 
were projected back to SU(3) elements. In test runs it turned out that varying the 
number of smearing iterations between 4 and 12 does not have a big impact on the 
length of the plateau in effective mass plots, see Figure [l|. Moreover, between the /3 
values analyzed we did not observe significant differences. 



3 Masses 

In this section we present our results on the masses. These were obtained from corre- 
lation functions Cr(p = 0,t), eq. at zero momentum in which case the exponential 
fall-off is given by the mass, Eh{p = 0) = Mh- For spatial correlation functions the 
exponential fall-off defines the screening mass. 

In order to obtain an estimate of the mass of the lowest state contributing to a 
given correlation function we first compared the effective mass plots, e.g. Figure II for 
various quark pair distances R in search of the optimal R value with regard to the onset 
and stability of the plateau. Subsequently, at the chosen R value fits over intervals 
[^min, N T — t min ] with varying t min (similarly for the spatial correlations at T ^ 0) were 
performed, again checking for stability of the mass value. Likewise, we symmetrized 
the correlation functions around midpoint, N T /2 and N a /2 respectively, and carried 
out fits including the full covariance matrix. Again, the minimum separation from the 
source was varied. The results quoted, mass values as well as errors, are obtained from 
the latter fits, selecting the fit interval by the best \ 2 value. The errors given include 
an estimate of the systematic error as suggested by a remaining dependence of the 
mass on R and t m [ n . Finally, also two-state fits were applied to correlation functions 
including data at next-to-optimal R values and smaller separations from the source in 
order to further check for consistency. 

The results for the ground state masses in the pseudoscalar and the vector channel 
at the temperatures investigated are summarized in Tables ^] to [| and are compared 
with the available zero temperature data in Figures |2| and |3]. In the figures we have 
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Figure 1: The effective mass M eS (t) as a function of t. The example shows the 
pseudoscalar channel at zero temperature at (3 = 6.2 and k = 0.141. The upper figure 
demonstrates the effect of varying the distance R between quark and antiquark while 
the number of smearing iterations was fixed to 4 in this example. In the lower plot, R 
is fixed at 4 and the number of smearing iterations is varied between and 12. 



plotted the so-called pole mass 



mu = 2 sinh 



M 



H 



(4) 
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Eq. (H) arises from a lattice meson action with a nearest neighbor symmetric difference 
discretization which appears to be favored by studies of lattice dispersion relations 
33, Using mu instead of the coefficient Mh of the exponential fall-off partly 



corrects for 0(am) lattice artefacts. We observe that my for the vectormeson at zero 
temperature, unlike My, is strikingly linear in the quark mass up to values in the 
vicinity of the charm quark mass. For small meson masses the difference between m# 
and Mh is of course negligible. 

Figure |^ summarizes the results at /3 — 6.0 which on the high temperature lattice 
corresponds to a temperature of 0.93 T c . The upper part of the figure shows my and 
m 2 P at small quark mass values while the lower part covers the entire quark mass 
range explored. The meson masses are plotted as function of m q where the quark 
mass is obtained from eq. ([!]) with k c (T = 0) determined from a linear fit in m q to 
the combined zero temperature literature data on mp p3| , |34|j up to quark mass values 



m q a < 0.08. As can be seen from the lowest lying line drawn in the upper part of 
Figure 0, the linear fit works well and deviations begin to emerge at m q a ~ 0.1. Our 
result for the critical kappa value is k c (T = 0) = 0.14542(2) which deviates marginally 
from the values quoted in |33],[34|, [3lJ because we have chosen a fit procedure slightly 
different from the one adopted there. Note however, that our value is well within 
the spread of the quoted k c values. Moreover, the upper part of Figure |] shows that 
the finite temperature pion (screening) mass retains a small but non-vanishing value 
at k c (T — 0). A fit to the mp data at quark mass values up to 0.08 of the form 
c + sm q , shown as the second lowest dotted line, leads to an intercept c = 0.006(3) at 
m q = 0. The slope s = 3.07(8) is a little larger than the value of s = 2.87(3) obtained 
from the equivalent fit to the T = data. Alternatively, one can perform a fit with a 
temperature dependent k c leading to a slightly different value, n c (T) = 0.14550(4). 

For the vectormeson masses, when plotted as a function of the zero temperature 
quark mass, we find that, as a whole, they apparently tend to be somewhat larger 
than the zero temperature values. In both cases, fits linear in the quark mass work 
well and we obtain my = 0.42(2) at T ~ 0.93 T c as opposed to 0.38(1) at T = in the 
chiral limit. This is shown as the appropriate lines in the upper part of Figure At 
least part of this difference could, however, be absorbed into the shift in k c mentioned 
above. 

In the lower part of Figure ^ my and m p are shown linearly over the entire quark 
mass range explored. At zero temperature, the vectormeson mass is linear in m q over 
the entire range. A fit ansatz including a term quadratic in the quark mass returns 
a value of 0.05(19) for its coefficient which is compatible with 0. The intercept is 
obtained as 0.381(6). Within errors this intercept is in agreement with the result of 
a linear fit, my = 0.379(2) + 2.588(5) m q which is shown in the figure. Note, that 
these numbers are in agreement with the results for the intercept obtained from the 
fit to the small quark mass data only. The mass of the pseudoscalar shows a behavior 
proportional to Jm q at small quark masses, as expected from chiral perturbation 
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16 3 x 32, /3 = 6.0, T = 



«1 


«2 


M P 


My 


0.1280 
0.1300 
0.1410 


0.1280 
0.1300 
0.1410 


1.259 (3) 
1.162 (3) 
0.559 (3) 


1.284 (4) 
1.191 (3) 
0.630 (4) 


24 3 x 8, /3 = 6.0, T = 0.93T C 


«i 




M P 


My 


0.1280 
0.1300 
0.1410 


0.1280 
0.1300 
0.1410 


1.262 (7) 
1.168 (6) 
0.557 (8) 


1.29 (1) 
1.197 (8) 
0.64 (1) 


32 3 x 8, (3 = 6.0, T = 0.93T C 






M P 


My 


0.1410 


0.14100 
0.14400 
0.14450 
0.14500 


0.559 (5) 
0.463 (6) 
0.443 (5) 
0.42 (1) 


0.64 (1) 
0.59 (1) 
0.57 (1) 
0.54 (1) 


0.1440 


0.14400 
0.14450 
0.14500 


0.323 (5) 
0.300 (1) 
0.27 (1) 


0.51 (2) 
0.50 (2) 
0.47 (2) 


0.1445 


0.14450 
0.14500 


0.273 (7) 
0.22 (1) 


0.46 (2) 
0.44 (3) 


0.1450 


0.14500 


0.20 (1) 


0.44 (2) 



Table 2: Meson masses Mh of the pseu- 
doscalar and vector meson at f3 = 6.0, 
T = oi2al6 3 x32 lattice and T = 0.93 T c 
on 24 3 x 8 and 32 3 x 8 lattices. 



24 3 x 48, (3 = 6.2, T = 





«2 


M P 


My 


0.1280 
0.1300 
0.1410 


0.1280 
0.1300 
0.1410 


1.124 (2) 
1.018 (2) 
0.333 (2) 


1.141 (2) 
1.040 (2) 
0.420 (4) 



24 3 x 8, (3 = 6.2, T = 1.23T C 


Ki 




M P 


My 


0.1280 
0.1300 
0.1410 


0.1280 
0.1300 
0.1410 


1.110 (4) 
0.995 (5) 
0.350 (4) 


1.110 (5) 
1.011 (4) 
0.410 (8) 



24 3 x 8, (3 = 6.2, T = 1.23T C 





^2 


M P 


My 


0.12800 


0.12800 
0.13600 
0.14151 
0.14232 
0.14280 


1.20 (2) 
1.04 (1) 
0.96 (1) 
0.96 (2) 
0.95 (2) 


1.23 (2) 
1.07 (1) 
1.01 (2) 
1.01 (2) 
1.00 (2) 


0.13600 


0.13600 
0.14151 
0.14232 
0.14280 


0.86 (1) 
0.78 (1) 
0.76 (2) 
0.76 (1) 


0.91 (2) 
0.85 (2) 
0.85 (3) 
0.85 (3) 


0.14151 


0.14151 
0.14232 
0.14280 


0.66 (2) 
0.65 (2) 
0.64 (2) 


0.76 (2) 
0.76 (2) 
0.76 (3) 


0.14232 


0.14232 
0.14280 


0.64 (2) 
0.64 (2) 


0.76 (2) 
0.75 (3) 


0.14280 


0.14280 


0.63 (1) 


0.75 (3) 



Table 3: Meson masses M H of the pseu- 
doscalar and vector meson at (3 = 6.2, 
T = on a 24 3 x48 lattice and T = 1.23 T c 
on a 24 3 x 8 lattice. 
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24 2 x 64 x 8, /3 = 6.4, T 


= 1.63T C 






M P 


My 


0.13000 


0.13000 


1.055 (5) 


1.085 (2) 


0.14000 


0.14000 


0.720 (4) 


0.756 (7) 


0.14030 


0.14030 


0.721 (4) 


0.757 (4) 


0.14060 


0.14060 


0.719 (3) 


0.756 (3) 


0.14090 


0.14090 


0.718 (3) 


0.755 (4) 



Table 4: Meson masses Mh of the pseu- 
doscalar and vector meson at /3 = 6.4, T = 
1.63 T c on a 24 2 x 64 x 8 lattice. 



theory, which turns into a linear one at larger m q values. Correspondingly, we chose a 
fit ansatz of the form 



m P = J bm q + dm 2 q (5) 

which returns b = 2.57(1) and d = 5.39(6), shown as the lowest line in the figure. 

Regarding the results at non-vanishing temperature, also at large quark masses 
we do not observe significant differences to the zero temperature data. However, the 
meson mass values are again plotted at the zero temperature quark masses, i.e. as a 
function of 1/k — 1/k c (T = 0) with our value of k c at T = 0. In a fit to the pseudoscalar 
mass data over the entire m q range, we therefore allow for a non-vanishing chiral limit 
and replace m q by m q — > m q + c in the ansatz eq. (|5|). This fit leads to the second 
lowest dotted curve in the lower part of Figure [I| with c = 0.005(2), b = 2.3(1) and 
d = 5.9(2). When fitting the vectormeson over all quark masses we observe a slightly 
better performance in terms of \ 2 °f a fit including terms quadratic in m q , my = 
0.413(6) + 2.45(5) m q + 0.13(5) m 2 ,, over a linear one, m v = 0.401(5) +2.59(1) m q . The 
quadratic term is not large though and the intercepts have their error bars touching 
each other. Note however that the intercept of the quadratic fit is in slightly better 
accord with the result of the linear fit to the small quark mass data. The corresponding 
line in Figure |^ is the quadratic function. 

Summarizing the results at T ~ 0.93 T c , we note that the pion (screening) mass 
exhibits only a very small value in the zero temperature chiral limit which could be 
absorbed by a slight shift in the critical kappa value. This might also be a finite 
size effect as the finite temperature lattice is somewhat bigger in spatial volume than 
the zero temperature ones. A non-vanishing temperature effect is perhaps seen in the 
vectormeson data which shows an enhancement of 7 to 10 % over the zero temperature 
data. It remains to be seen, however, whether this finding persists in the continuum 
limit. On the other hand, our results in the Wilson discretization are in good agreement 
with observations made in earlier studies using staggered fermions 1T3 . 
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Figure 2: Meson masses at T ~ 0.93 T c compared with the zero temperature results 
as a function of the quark mass. The zero-temperature data is partially taken from 
the literature |33], [34], |39|. The upper plot features the vectormeson mass and the 
pseudoscalar mass squared in the vicinity of the chiral limit while in the lower part 
the data for the entire quark mass range are shown. Here both masses are plotted 
linearly. The lines are fit results explained in the text. 
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Figure 3: Meson masses at T ~ 1.23T C and 1.63T C compared with zero temperature 



data partially taken from the literature |33l, j34|, |39 . 



The results for the meson masses at (3 = 6.2 and 6.4 corresponding to the two tem- 
peratures above T c which are investigated here, T ~ 1.23 T c and 1.63 T c respectively, 
are summarized in Figure || Again, we show my and mpasa function of the quark 
mass m q with k c (T = 0) determined as explained above from fits to the combined 
zero temperature pion mass results [|34],|35| at small quark mass. The resulting critical 
k values are k c = 0.14313(2) and 0.14141(3) at 6.2 and 6.4, respectively. 
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The zero temperature data shows the same behavior as at (3 = 6.0, including the 
linearity in m q of the vector meson mass up to the largest quark masses explored. At 
(3 = 6.2 we fitted my at T = over the entire quark mass range both, quadratically 
in m q , leading to my = 0.283(2) + 2.61(2) m q + 0.10(4) m q as well as linearly, my = 
0.279(1) + 2.65(5) m q . Again, the quadratic term was found to be small and the 
intercepts are stable. Using only data close to the chiral limit leads to an extrapolated 
value of m p = 0.267(8), so that the deviation from the fits to the full data set is a 
little larger than at (3 = 6.0. The dotted lines given in Figure || are the quadratic fits. 
Note, that the data points labelled by "non-deg" indicate meson masses obtained 
from propagators with two quarks of non-degenerate masses. They are plotted as a 
function of an "effective" quark mass m q = (m q> i + m q ^)j2 and fit very well on the 
quadratic curves at T = 0. The pseudoscalar data was again fitted with the ansatz of 
eq. (|5|) which gives b = 1.846(8) and d = 6.22(3) in this case. This is shown again as 
the lowest line in the figure. 

At (3 = 6.4 only data at small quark masses were available |$3J and we have 
carried out a linear fit to the vector meson, with the result my = 0.19(1) + 3.23(3) m q . 
Although k c (T = 0) was obtained from fits to m 2 P in the vicinity of the chiral limit, in 
Figure ^| we only show, as lines to guide the eye, the results of fits with eq. @ under 
the condition that mp <my for large m q . 

In contrast to the situation at T < T c , at both temperatures above the transition 
point the screening masses are markedly different from the zero temperature results, 
especially at small quark mass values. The curves interpolating the non-zero tem- 
perature data are polynomials of third degree for the vectormeson at T = 1.23T C and 
of second degree elsewhere. We have chosen this ansatz because above T c no par- 
ticular functional form is known so far. The curves are the results of fits to data 
for degenerate quark masses and are mainly meant to guide the eye. The intercepts 
were obtained from these fits as m 7r (1.23T c ) = 0.63(2), m p (1.23T c ) = 0.74(2) and 
77^(1.637^) = 0.735(5), m p (1.63T c ) = 0.774(8) at the two temperatures respectively. 
They are discussed in the following. 

At high temperatures one expects that the plasma consists of only weakly inter- 
acting quarks and gluons. Correspondingly, the correlation functions with mesonic 
quantum numbers should be described to first approximation by the free propagation 
of a quark-antiquark pair. In this case the exponential fall-off of the correlation func- 
tions is not dominated at large distances by a single mass originating from an isolated 
pole in the spectral function but rather by all possible energies of quark and antiquark 
with opposite "momenta" of equal size adding to zero total momentum to which the 
correlation functions are projected on. Correspondingly, effective masses taken from 
point to point correlations (R — in eq. (§)) exhibit a rather marked curvature as 
function of the separation between the meson operators. Tuning R away from leads 
to a suppression of higher quark momenta in the momentum sums. This way, it is pos- 
sible to obtain plateau-like behavior in effective mass plots also in the free quark case. 
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The lowest value of quark "momentum" in the temporal direction, i.e. the smallest 
Matsubara frequency is ttT due to the antiperiodic boundary condition for fermions in 
the t direction. This leads to the lowest "energy" contributing to spatial correlation 
functions of mesonic operators of 



m H = 2^m\ + (ttT) 2 (6) 

in the continuum. On a finite lattice, eq. ([]) is subject to finite volume and finite 
lattice spacing corrections, see e.g. J7|. In order to compare our simulation results at 
T > T c with the free quark case, we have computed the free quark propagator on a 
finite lattice analytically and performed the sum over all lattice "momenta" p of the 
quarks numerically. For definiteness, we then plot the effective mass for point to point 
correlations at a separation z = N a /A. This results in the solid curve shown in the 
upper part of Figure [|. 

As in the free case, in our results from the Monte Carlo simulations we are also 
able to identify plateaux in the effective mass plots when the separation R, eq. (0), 
is tuned. The vectormeson and, to a somewhat less degree, the pseudoscalar masses 
approach the free quark limit as is seen from Figure []. The remaining differences are 
becoming smaller with rising temperature. Moreover, the fact that the masses for non- 
degenerate quark combinations, plotted at m q = (m qt i + m q ^)/2, are above the ones 

for degenerate quarks would fit into this picture as 2^rn 2 + (vrT) 2 < J rn 2 ql + (ttT) 2 + 



m q2 + (jtT) 2 . At large quark masses the Monte Carlo data exceed the free quark 

curve. This could be explained by noting that in the lattice version of eq. @ we have 
used the bare quark mass. At the temperatures investigated one should presumably 
compare with an effective quark mass which also accounts for a thermal contribution 
||. In addition, it has been argued in the context of dimensional reduction that the 
(confining) potential of the reduced three-dimensional theory leads to modifications 
of eq. (H) for screening masses pi . These modifications will be positive and of order 



^/(Jspat, where <7 spat is the string tension of spatial Wilson loops, and can be estimated 



to amount to an 0(10%) effect ^3]. 



Finally, in Figure [| our results are compared with data obtained in the staggered 
fermion discretization, both below T c |12|] and above [|, |7|. Below T c the 7r and p 
screening masses are shown in the chiral limit. The agreement between the two dis- 
cretizations is evident. Above T c the staggered data was computed at a small bare 
quark mass of m q a = 0.02 while the Wilson results have been extrapolated to the chiral 
limit. From Figure 0, however, recall that above T c the screening masses are practi- 
cally independent of m q at small quark masses. Both data sets have been rescaled by 
the appropriate ratio of 2ir, eq. (|D, to its finite lattice size corrected value so that the 
data can immediately be compared with the continuum expectation for free quarks. 
Figure || shows that above T c the p screening mass rapidly approaches 2ir. This holds 
in both discretizations. The pion however behaves very differently for the two lattice 
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actions. In the Wilson case, at 1.63 T c the pion screening mass is almost degenerate 
with the p and close to 2tt, with a it to p ratio of 0.95(1). The same approximate 
degeneracy of it and p screening masses has been observed in the only other quenched 
analysis with Wilson quarks available so far JTTj . Here a ratio of 0.955(7) has been ob- 
tained at a temperature of about 1.5 T c . Quite contrary, the staggered pion screening 
mass is much smaller than the p, m n /m p = 0.75(2) at 1.8T C , and even at 5T C reaches 
only about 75 % of two times the lowest Matsubara frequency. There is no immediate 
explanation at our hands at least. Possible reasons include the different symmetries 
of the lattice actions at non-vanishing lattice spacing although the lattice spacings are 
below 0.1 fm, the different aspect ratios although it is hard to see why this should 
affect the 7r to p ratio, and finite volume effects although those would be expected to 
go the opposite way. Thus, it is hoped that future systematic studies help to resolve 
this discrepancy between the two lattice fermion formulations. 




Figure 4: Comparison with the results of quenched staggered simulations below (12 
and above T c The data have been rescaled by the appropriate ratio of 2n, eq. (H), 
to its finite lattice size corrected value so that the data can immediately be compared 
with the continuum expectation for free quarks shown as the horizontal line. 



4 Dispersion relation 

At non-zero temperature, Lorentz invariance is broken because the temporal direction 
is distinguished as the direction of the four- velocity of the heat bath. As a consequence, 
unlike the zero temperature case where it depends on the Lorentz invariant scalar p 2 , 
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in this case the spectral density will depend on temporal and spatial components po 
and p separately. At temperatures below the confinement-deconfinement transition 
the spectrum will still consist of particle excitations, yet, their dispersion relations 
might be more complicated and reflect the breaking of Lorentzian invariance. The 
spectral density will be of the form 

p(p ,p) = 27re(po)5(pl-cu 2 (p,T)) (7) 

where 

uj 2 (p,T) = m 2 +p 2 + Tl(p,T) (8) 

contains the temperature dependent vacuum polarization tensor U(p,T). As a simple 
example, assume that the temperature effects can be absorbed into a temperature de- 
pendent mass m(T) and a coefficient A(T) which might also be temperature dependent 
and different from 1, 

uj 2 (p, T) ~ m 2 (T) + A 2 (T)p 2 (9) 

Such an approximation might hold at least at small temperatures. In this case, at zero 
momentum the temporal correlator will decay with the so-called pole mass m(T) 

C(p = 0,t) ~ exp(-m(T)t) (10) 

whereas the spatial correlation function has an exponential fall-off 

C0=O,z) ~exp(-m sc (T)z) (11) 

determined by the screening mass m sc (T) = m(T)/A(T) which differs from the pole 
mass if A(T) ^ 1. At non- vanishing "momentum", p ^ 0, the exponential decrease of 
the spatial correlator is described by u sc , 

C(p,z) ~ exp(-u sc z) (12) 

where in this particularly simple example u sc is given as 

U) 2 

"sc=Pl + ^+™ 2 c (13) 

Comparing different projections to Matsubara frequencies u n = 2irTn and to spatial 
momenta p± might thus reveal a dispersion relation different from the zero temperature 
one and, moreover, indicate a coefficient A different from 1 and correspondingly a 
difference between pole and screening mass. 

The formulae given above are the continuum dispersion relations. They are modi- 
fied on the lattice. Zero temperature studies have shown that a lattice disper- 
sion of the form 

sinh2 (f)=E- 2 (f)+^ 2 (f) (") 
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arising from an effective boson action with a nearest-neighbor kinetic term, is best 
capable to describe lattice data at non-zero momentum. 

We investigated the lattice dispersion relation of mesons below T c , T ~ 0.93 T c , 
for different quark masses and compared it with T = data. In Figure |^ we show 
our results obtained in the pseudoscalar channel at k — 0.141. We have plotted M^a 
obtained by subtracting the momentum contribution sin 2 (pfca/2) from both sides of 
eq. (|i"4]) where E is the fitted coefficient of the exponential fall-off in z of the spatial 
correlation function. Results at other k values investigated are very similar. 
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Figure 5: The lattice dispersion relation shown for the pseudoscalar channel at 
k = 0.141. In the figure, "dir" indicates whether spatial or temporal components 
of {pxiPyiPt) were chosen to be different from 0. The horizontal lines denote the dis- 
persion relation eq. (F3J). Note that the T = data has been shifted upwards, right 
scale. 



In Figure §, a pk independence of the data indicates that eq. fll4|) is the correct 
dispersion relation. Indeed, it seems to be favored by the zero temperature data. As 
in ||38|| , using a different dispersion relation led to data points rising or falling with p^. 



Regarding the data at T = 0.93 T c , again we observe that the zero-temperature 
dispersion, eq. (0), is describing the data for non- vanishing spatial momentum com- 
ponents properly. Note that the figure contains data from two different lattice sizes, 
leading to different values for the spatial momenta. Moreover, the single data point for 
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the lowest bosonic Matsubara frequency uj\ = 2irT is also lying on the horizontal line, 
suggesting that the coefficient A of eq. fllBp is not too far from unity. Of course, this 
conclusion is tied to the applicability of eq. fll3|) , moreover, the statistical significance 
is certainly not overwhelming. Nevertheless, this result might be taken as further 
support for the difference between pole and screening masses not being too large at 



T ~ 0.93 T c , see also |7 



In addition to the investigation below T c , we also computed the dispersion relation 
in mesonic channels in the deconfined phase. The results are shown and compared 
with the free quark case in Figure |^. The former have been determined by fits in 
the same way as the masses while for the latter, to be definite we use the effective 
energies at N a /4, obtained from sink operators with R = 4. Both, Monte Carlo as 
well as free quark results are somewhat dependent on the quark-antiquark distance. 
As can be read off the figure, at 1.23 T c the data for spatial momenta show systematic 
deviations from the preferred zero-temperature dispersion relation. These differences 
have almost disappeared at 1.63 T c and the data are becoming consistent with the zero- 
temperature dispersion relation which also happens to describe the free quark case. 
For the lowest temporal mesonic "momentum", p t = 2ttT, at 1.23 T c pseudoscalar and 
vector "energy" both are systematically lower than the corresponding results for a 
spatial momentum of exactly the same value showing that Lorentzian symmetry is 
disturbed. At 1.63 T c the data follow the free quark behavior, in particular insofar in 
the vector channel the "energy" at p t = 2ttT is much smaller than the "energy" at the 
same spatial momentum^. 



5 Wave function 

Bethe-Salpeter amplitudes provide information about the probability of finding a pre- 
arranged configuration of quarks inside a hadron. Their general definition is given 
by 

<S>(R) = {0\O(R)\H(f)) (15) 

where, for \H) being a meson state, the operator O(R) annihilates a quark-antiquark 
pair with the appropriate quantum numbers separated by R. The choice of the 
operator is not unique. On the lattice, Bethe-Salpeter amplitudes have been stud- 
ied in Coulomb and in Landau gauge or by using various gauge-invariant defini- 
tions []3~6"|, |40l-f43l . The various methods treat the gluon flux tube connecting quark 
and antiquark in a different way so that it may not be surprising that different defini- 
tions of the amplitudes have lead to different results (see however Moreover, an 

a Note that this, at first sight strange behavior can easily be reproduced in the continuum if only 
the lowest quark Matsubara frequency would contribute to the correlation function. Of course, in 
this case pseudoscalar and vector channel would lead to the same result. Apparently, the interplay 
between higher quark momenta and the spin structure is important. 
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Figure 6: The lattice dispersion relation at two temperatures above T c and in the free 
quark case. Shown are both the results for the pseudoscalar and the vector channel. 
The bare quark mass is roughly equal to m s . Again, "dir" indicates whether spatial or 
temporal components of {p x ,Py,Pt) were chosen to be different from 0. The horizontal 



lines denote the dispersion relation eq. ([14]) . 
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immediate connection between the size of a hadron as extracted from the amplitudes 
and the measured radii has not been demonstrated. Nevertheless, the Bethe-Salpeter 
amplitudes do provide some qualitative insight into the quark distribution inside a 
hadron, and in particular they should be able to capture differences in the structure 



of hadronic excitations at zero and high temperature if these are present [[19 



As a by-product of the attempt to obtain a better projection onto the ground 
state, in this paper we use a gauge-invariant definition of 0(R) as given in eq. (fj) f3"6fl . 
We analyze S-wave states only as this operator is invariant under 90° rotations. The 
Bethe-Salpeter amplitudes are determined from simultaneous two-state fits to all the 
correlation functions Cr at all R values at which they were calculated. In these fits 
the ground state mass was required to be the same at all R whereas the mass of the 
second state as well as the amplitudes for both ground and excited state were allowed 
to be R dependent. Again, the fit interval was varied and we checked for stability of 
the fit results under this variation as well as for agreement of the ground state mass 
with the numbers obtained as explained in section |^. 

The procedure proved stable enough to extract the Bethe-Salpeter amplitude for 
the ground state as well as to obtain an estimate for the excited state. An example of 
the results is shown in Figure 0. Here, the wave functions have been normalized to 1 
at R = 0. Note that the amplitude of the excited state is vanishing at about 1/2 fm. 
This corresponds with the observation that at large R effective masses approach the 
plateaux in effective mass plots from below, see Figure [l|. The results for the excited 
states are more sensitive to variations of the fit interval since at small separations from 
the source more than one excited state can contribute whereas at large separations 
the excited state dies out. We will thus concentrate on the ground state results in the 
following. 

The ground state wave functions Q are fitted with slightly modified hydrogen S-wave 



functions 

<$>(R) oc exp(-(R/a ) u ) (16) 

with the free parameters ao and v. In non-relativistic potential models the power of 
R is obtained as v = 1 for a purely Coulombic and v = 3/2 for a linear potential. The 
wave functions can be used to compute an estimate of the size of the meson , 



= / i'R if ) W)/ / d'S. *>(R> = ffirSuy (17) 



The results of these fits to the pseudoscalar as well as vector meson ground state 
are summarized in Table |^. We first want to stress that the measured zero temperature 
wave functions, i.e. the parameter ao and the radius respectively, scale with (3 as the 

b The excited states could be well fitted by the ansatz oc (1 — R/bi) exp(— (R/aiY) which 

accounts for the node in the wave function. 
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Figure 7: Wave functions for the ground (GS) and excited state (ExS) of the pseu- 
doscalar at T = and T = 0.93T C , for m q K,m s 



lattice spacing when the (physical) quark mass is fixed. This is lending further support 
that the investigated observable is physical. 

At T = 0, the obtained values for the parameter v lie between the values for a purely 
Coulombic and a purely linear potential. The values for v are slightly increasing with 
the quark mass and they are larger for the vector meson than for the pseudoscalar. 
Moreover, the values for and for the radii are also increasing with the quark mass 
and again are larger for the vector meson than for the pseudoscalar. Both effects are 
consistent with the experimental observation that the p meson is larger than the pion 
and thus probes the linear rising part of the potential at larger distances. Compared 



to the (Coulomb gauge) results of |42| at the same pseudoscalar to vector meson 
mass ratio, we observe perfect agreement for the parameter v of the pseudoscalar 
ground state whereas our v value for the vector meson is slightly smaller. The second 



moments of the meson wave functions, eq. flT7D, are about 10 % larger than in [12 
Still, as reported there, we obtain radii which are about a factor two smaller than 
sizes obtained from form factor measurements. Whether this discrepancy is due to the 
quenched approximation i.e. to not taking proper account of a pion cloud still has to 
be speculated at the present stage. 

At T = 0.93T C the pseudoscalar as well as vector meson (spatial) wave functions 
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Table 5: Results for the coefficients a and v of the ground state wave function, 
eq. (0), as well as for the mean radius squared, eq. (|T7|), for the pseudoscalar (P) and 
vector meson (V). 



agree with the (temporal) zero temperature ones at quark masses where we have data 
for both. This also holds for the excited state, see Figure 0. At this temperature we 
also have data at small quark masses so that we can attempt a chiral extrapolation. 
A fit to the squared radii which is linear in m 2 PS leads to values of ((r/a) 2 ) = 7.6(6) 
and ((r/a) 2 ) = 9.6(12), respectively, which are in agreement within errors with the 



zero-temperature numbers reported in [42 



At the two temperatures above T c it is observed that the spatial wave functions in 
lattice units practically remain unchanged between the two (3 values. This is distinc- 
tively different from the zero temperature situation. 



In more detail, fits to the wave functions with eq. ( flq) return the values given in 
Table ||. The results for the parameter clq are almost independent of (3. The parameter 
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v takes on values which are larger than at zero temperature. Close to the chiral limit, 
v changes from about 1.4 at 1.23T C to 1.6 at 1.63T C for the pion while for the p the 
change with temperature is smaller, u ~ 1.8 and 1.9 respectively. The second moment 
of the wave function, (r 2 ), in lattice units stays constant within errors for the p while 
for the pion with rising temperature we observe a decrease of the central values, with 
error bars just overlapping. As the quark mass is increased the size of the state is 
slightly shrinking. 

Translated to physical distances a shrinking of the states' size with temperature is 
clearly oberved, Figure [|. At least qualitatively, this observation is in accord with the 
anticipated behavior [p0|,|2~4j| of 



$(R) oc exp(- Ja s UT)7rTR 3 / 2 ) 



where cr spat (T) is the temperature dependent string tension of spatial Wilson loops, 
\J<J spaX {T) = cg 2 (T)T. Given the small variation of the effective coupling constant 
g(T) between T = 1.23 and 1.63T C , (r 2 ) should thus shrink proportional to 1/T 2 to 
first approximation. Indeed, this seems to hold for the vector meson whereas the 
variation with temperature of the pion radius is stronger. Moreover, the power of the 
R dependence is not quite as anticipated. A more detailed investigation along the 
lines of p0| is planned for the future. 
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6 Conclusion 



In this paper we have presented results of a lattice investigation of meson correlation 
functions built from improved Wilson quarks on quenched gauge field configurations. 
Below the confinement-deconfinement transition temperature we do not observe any 
significant difference between finite temperature screening masses and zero tempera- 
ture masses. This agrees with earlier findings obtained in the staggered discretization. 
Likewise, the spatial wavefunctions at non-vanishing temperature agree with wave 
functions at T = 0. Moreover, the dispersion relation does not seem to depend on 
whether the "momenta" are chosen to point in the spatial or temporal direction. This 
can be viewed as a hint towards the similarity of screening and pole masses as has also 
been reported in an earlier study on anisotropic lattices. Above the critical temper- 
ature, we have computed meson screening masses which in both cases, pseudoscalar 
as well as vector channel, are close to twice the lowest quark Matsubara frequency. 
Moreover, at 1.63 T c the "momentum" dependencies of the correlation functions are 
in agreement with the free quark case while at 1.23 T c the situation certainly is more 
complex. Finally, the spatial wave functions, at least qualitatively, behave in accord 
with the expectation that the temperature dependence of the spatial string tension 
leads to a narrowing of the wave functions with rising temperature. We thus come to 
conclude that above T c the spatial correlation functions are consistent with a leading 
two quark contribution. This still leaves room for infrared non-perturbative effects. 
Note in this respect that HTL-resummed perturbation theory does not show significant 
changes in (temporal) correlation functions although the underlying spectral function 
is quite different from the free case | 44| . It is our hope that more refined analysis 



techniques along the lines set in [45f and perhaps more systematic studies of infrared 
properties on large lattices extend the knowledge about hadronic excitations in the 
plasma phase. 
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